function [nmse, rmse] = compute_NMSE_RMSE(error, vali)
rmse = rms(error,2);
nmse = rmse./rms(vali,1)';
end